Analysing polynucleotide sequences ^ 

, ^^This is a CIP of Serial No. 07/573,317 f iled .J 
September 28, 1990. . -- " 

1. INTRODUCTION 

Three methods dominate molecular analysis of 
nucleic acid sequences: gel electrophoresis of 
restriction fragments, molecular hybridisation, and the 
rapid DNA sequencing methods. These three methods have 
a very wide range of applications in biology, both in 
basic studies, and in the applied areas of the subject 
such as medicine and agriculture. Some idea of the 
scale on which the methods are now used is given by the 
rate of accumulation of DNA sequences which is now 
well over one million base pairs a year. However, 
powerful as they are, they have their limitations. The 
restriction fragment and hybridisation methods give a 
coarse analysis of an extensive region, but are rapid; 
sequence analysis gives the ultimate resolution, but it 
is slow, analysing only a short stretch at a time. ... 
There is a need for methods which are faster than the 
present methods, and in particular for methods which 
cover a large amount of sequence in each analysis. 

This invention provides a new approach which 
produces both a fingerprint and a partial or complete 
sequence in a single analysis, and may be used directly 
with complex DNAs and populations of RNA without the 
need for cloning. 

In one aspect the invention provides 
apparatus for analysing a polynucleotide sequence, • 
comprising a support and attached to a surface therof 
an array of the whole or a chosen part of a complete' 
set of oligo nucleotides of chosen lengths, the 
different oligonucleotides occupying separate cells of 
the array and being capable of ta k i ng- 'part in 
hybridisation reactions. For studying differences 
between polynucleotide sequences, the invention- 
provides in another aspect apparatus comprising a 
support and attached to a surface therof an array of 
the whole or a chosen part of a complete set of 
oligonucleotides of chosen lengths comprising the 
polynucleotide sequences, the different 



oligonucleotides occupying separate cells of the array 
and being capable of taking part in hybridisation 
reactions . 

In another aspect, the invention provides a 
method of analysing a polynucleotide sequence, by the 
use of a support to the surface of which is attached an 
array of the whole or a chosen part of a complete set 
of oligo nucleotides of chosen- lengths, the different 
oligonucleotides occupying separate cells of the array , 
which method comprise-s labelling the polynucleotide 
sequence or fragments thereof to form labelled 
material, applying the labelled ma ter i a 1 under 
hybridisation conditions to the array , a'nd'observ i ng 
the location of the label on the surface associated 
with particular members of the set of oligonucleotides. 

The idea -of the invention is thus to provide 
a structured array -of the whole or a chosen part of a 
complete set of oligonucleotides of one or several 
chosen lengths. The array,. which may be laid out on a 
supporting film or glass plate, forms the target for a 
hybridisation reaction. The chosen conditions of 
hybridisation and the length of the oligonucleotides 
must at all events be sufficient for the available 
equipment to be able to discriminate between exactly 
matched and mismatched oligonucleotides. In the 
hybridisation reaction, the array is explored by a 
labelled probe, which may comprise oligomers of the - 
chosen length or longer polynucleotide sequences- or 
fragments, and whose nature depends on the particular 
application. For example, the probe m : ay comprise 
labelled sequences amplified from genomic DMA by the 
polymerase chain reaction, or a mRNA population, or a 
complete set of oligonucleotides from a complex 
sequence such as an entire genome. The end result is a 
set of filled cells correspond i ng.. to the 
oligonucleotides present in the analysed sequence, and 
a set of "empty" sites corresponding to the*, sequences 



which are absent in the analysed sequence. The pattern 
produces a fingerprint representing ail of the sequence 
analysed. In addition, it is possible to assemble aost 
or ail of the sequence analysed if an oligonucleotide 
5 length is chosen such that most or all oligonucleotide 
sequences occur only once. 

The number the length and the sequences of the 
oligonucleotides present in the array "lookup table" 
also depend on the application. The array may include 
10 all possible oligonucleotides of the chosen length, as 
would be required if there was no sequence information 
on the sequence to be analysed. In this case, the 
preferred length of oligonucleotide used depends on the 
length of the sequence to be analysed, and is such that 
13 there is likely to be only one copy of any_ particular 
oligomer in the sequence to be analysed. Such arrays 
are large. If there is any information available on 
the sequence to be analysed, the array may be a 
selected subset. For the analysis of a sequence which 
2G is known, the size of the array is of the sane order as 
length of the sequence, and for rany applications, such 
as the analysis of a gene for mutations, it can be 
ruite small. These factors are discussed in detail in 
what follows. 
25 2. OLIGONUCLEOTIDES A3 SEQUENCE ??.03ES 

Oligonucleotides form base paired duplexes with 
oligonucleotides which have the complementary base 
sequence. The stability 'O f the duplex is dependent on 
:he length of the oligonucleotides and on base 
3 0 oo ^position. Effects of base exposition on duplex 
stability can be greatly reduce: by the presence of 
high concentrations of quaternary or tertiary amines. 
However- there is a strong effe:t of ~iz matches in tne 
:ligc nucleotides duplex on tne :'r.er:ai stability cf -he 
' 35 r.vbrid and it is this which ma>es the technique of 
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hybridisa tion with oligonucleotides such a powerful 
method for the analysis of mutations, and for the 
selection of specific sequences for amplification by 
D it A polymerase chain reaction. The position of the 
mismatch affects the degree of des tabilisa tion . 
Mismatches in the centre of the duplex may cause a 
lowering of the Tm by 10°C compared with 1°C for a 
terminal mismatch. There is then a range of 
discriminating power depending on the position of 
mismatch, which has implications for the method 
described here. There are ways of improving the 
discriminating power, for example by carrying out 
hybridisation close to the Tm of the duplex to reduce 
the rate of formation of mismatched duplexes, and by 

1 " increasing the length of oligonucleotide beyond what is 

required for- unique representation. A way of doing 

this systematically is discussed. 

3. ANALYSIS OF A PREDETERMINED SEQUENCE 

One of the most powerful uses of oligonucleotide 
20 probes has been in the detection of single base changes 
in husar. genes. The first example was the detection of 
the single base change in the betaglocin gene which 
leads to sickle ceil disease. There is a need to 
extend this approach to genes in which there may be a 

2 5 number of different mutations leading to the same 

pheno:y:e , for example the DMD gene arc the HPRT gene, 
and to find" an efficient way of scanning the human 
genome for mutations in regions which have been shown 
by linkage analysis to contain a disease' locus for 

30 example Huntington's disease and Cystic Fibrosis. Any 
known sequence can be represented completely as a set 
cf overlapping oligonucleotides. The size of the set 
is N s + 1 N , where N is the lengtr. c f the sequence 
and s is the length cf an oligomer. A gene of 1 k b for 

35 example, may be divided into an overlapping "set of 
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around one thousand oligonucleotides of any chosen 
length. An array constructed with each of tr.ese 
oligonucleotides in a separate cell can be used, as a 
DU ltiple hybridisation probe to examine the homologous 
sequence in any context, a single-cooy gene in the hu,an 
genome or a messenger RNA among a mixed RNA peculation, 
for example. The length s may be chosen such that 
tr-r» is only a small probability that any oligomer in 
the sequence is represented elsewhere in the sequence 
to be analysed. This can be estimated from the 
exoression given in the.se-ction discussing statistics 
below. For a less complete analysis it would ce 
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.bie tojeducs the size of the array e.g. by a 



factor of up to 5 by representing the sequence in a 
15 partly or non-overlapping set. The advantage of using 
a comoletely overlapping set is that it provides a more 
precise location of any sequence difference, as the 
mismatch will scan in s consecutive oligonucleotides. 
If. ANALYSIS OF AN UNDETERMINED SEQUENCE 
20 The genomes of ail free living organisms are 

larger than a million base pairs and none has yet been 
seouenced completely. Restriction site matting reveals 
only a small tart of the sequence, and can detect only 
a small portitr. of mutations when used to compare two 
25 genomes. Mere efficient methods for analysing complex 
sequences are needed to bring the full power :f 
molecular genetics to bear on the many biological 
problems for vr.ich there is no direct access to the 
gene cr genes involved. In many cases." the full 
30 seouer.ee of tr.e nucleic acids need not be de t e r mine d ; 

the important sequences are those whic.n differ between 
two nucleic a:ids. To give three examples: tr.e DNA 
sequences vhi:n are different between a wild type 
organism ar.i :.-.e which carries a mutant car. lead the 
way to isolation of the relevant gene; similarly, the 
sequence differences between a cancer cell and its 
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normal counterpart car. reveal the cause of transfor- 
mation; and the R N A sequences which differ between :v: 
ceil types point to the functions which distinguish 
them. These problems can be opened to molecular 
analysis by a set hod which identifies sequence 
differences. Using the approach outlined here, such 
differences can be revealed by hybridising the two 
nucleic acids, for example the genomic DNA of the two 
genotypes, or the nKNA populations of two cell types :c 
an array of oligonucleotides which represent all 
possible sequences. Positions in the array which are 
occupied by one sequence but not by the other show 
differences in two sequences. This gives the sequence 
information needed tc synthesise probes which can then 
15 be used to isolate clones of the sequence involved. 

4.1 ASSEMBLING. THE SEQUENCE INFORMATION 
Sequences can be reconstructed by examining the 

result of hybridisation to an array. Any oligonuclec- 
^ tide of length s frcz within a long sequence, overlaps 

^W$$ 20 with two owners over a length s-K Starting from each 

positive Oligonucleotide , the array may be examined : 2 : 
j=i the four oligonucleotides to the left and the four t c 

tfl the right that can overlap with a one base displaces r. * 

'If only one of these four oligonucleotides is found tc 
25 be positive to the rignt, then the overlap and the • 
additional. base to the right determine s bases in the 
unknown sequence. The process is repeated in both 
directions, seeking unique matches with other positive 
oligonucleotides in t.te array. Each unique natch = 
30 a base to the reconstructed sequence. 

4 . 2 SOME STATISTICS 
Any sequence c f length N can be broken down t c 1 

set of ~ M overlapping sequences s base pairs in 
0^ length. (For double stranded nucleic acids, the 

35 sequence complexity z: ^ sequence of U base pairs is 
2N , because the two strands have different sequences. 
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but for the present purpose, this factor of two is not 

significant). For oligonucleotides of length s, there 
. s 

are 4 airrerent sequence cor b i nations. How big should 
s be to ensure that most oligonucleotides will be 
represented only once in the sequence to be analysed, 
of complexity N? For a random sequence the expected 
number of s-mers which will be present in more than one 
co p y i s 

/i>i *4'(1 + A)) 

where 

A = (vY-j-f l)/4' 



For practical reasons it is also useful to know 
how many sequences are related to any given s-mer by a 
single base change. Each position can be substituted 
by one of three bases, there are therefore 3s sequences 
related to an individual s-:er by a single base change, 
and the probability that any s-mer in a sequence of N 
bases is related to any other s-mer in that sequence 
allowing one substitution is 3s x N/M S . The relative 
signals of. matched and mismatched sequences will then 
depend on hew good the hybridisation conditions are in 
distinguishing a perfect xa;:h from one which 'differ by 
a single base. (If h~ is an order of magnitude greater 
than N, there should only be a few, 3s/10, related to 
any oligonucleotide by one base change.) The 
indications are that the yield of hybrid froa the mis- 
matched sequence is a fraction of that formed by the 
perfect duplex. 

For what follows, it i z assumed that conditions 
can be found which allow :lir:nu:iectices which have 
complements in the probe to be d i s t i ngu i s n e d from those 
v h i c h do not. 

4*3 ARRAY FORMAT, CONSTRUCTION" A N D SIZE 

To form an idea of the scale of tne arrays needed 
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to analyse sequences of different co-plexity it is 
convenient to think of the array as a square matrix. 
All sequences of a given length can be represented just 
once in a matrix constructed by drawing four rows 
representing the four bases, followed by four similar 
columns. This produces a 4 x 4 matrix in which each of 
the 16 squares represents one of the 16 doublets. "our 
similar matrices, but one quarter the size, are then 
drawn within each of the original squares. This 
produces a 16 x 16 matrix containing all 256 tetra- 
nucleotide sequences. Repeating this process produces 
a matrix of any chosen depth, s, with a nuaber of cells 
equal to 4 s . As discussed above, the choice of s is c: 
great importance, as it determines the complexity of 
the sequence representation. As discussed below, s 
also determines the size of the matrix constructed, 
which must be very big for complex genomes. Finally, 
the length of the oligonucleotides determines the 
, hybridisation conditions and their discriminating power 
as hybridisation probes. 
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The table shows the expected scale cf the arrays 
needed to perform the first analysis of a few genomes. 
The examples were chosen because they are genomes which 
have either been sequenced by conventional procedures - 
5 the cosmic scale - , are in the process of being 

sequenced - the E. coli scale or for .which there has be; 

considerable discussion of the magnitude of the problem 
- the human scale. the table shows that the expected 
scale of the matrix approach is only a small fraction 
10 of the conventional approach. This is readily seen in 
the area cf X-ray film that would be consumed. It is 
also evident that the time taken for the analysis would 
be only a snail fraction of that needed for gel 
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,n methods. The "Genomes" column shows the length of 

nj . 15 random sequence which would fill about'5% of cells in 

the matrix. This has been determined to be the optimum 
condition for the first step in the sequencing strategy 
OR discussed below. At this size, a high proportion of 

the positive signals would rep re sent single occurrences 
20 of each oligomer, the conditions needed tc compare two 
fj genomes for sequence differences. 

MB 5. REFINEMENT Or AN INCOMPLETE SEQUENCE 

Reconstruction of a complex sequence produces a 
result in which the reconstructed secuence is 
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by these repeats can be removed by extending the 
sequence to longer oligomers. In principle, those 
sequences shown to be repeated by the first analysis, 
using an array representation of all possible 
5 oligomers, could be r e s yn the s is ed with an extension at 
each end. For each repeated oligomer, there would :e 4 
x 4 = 1 6 oligomers in the new matrix. The 
hybridisation analysis would now be repeated until tne 
sequence was complete. In practice, because the 
10 results of a positive signal in the hybridisation = a y 

be ambiguous, it may be better to adopt a refinement of 
the first result by extending ail sequences which tit 
O not give a clear negative result in the first analysis. 

4f An advantage of this approach is that extending 



m 



Id sequence brings ciisratches whicn are close to tne er.zz 
Ln in the shorter oligomer, closer to the centre in tne 

0 1 extended oligomer, increasing the discriminatory pover 

of duplex formation. 

5.1 A HYPOTHETICAL ANALYSIS OF THE SEQUENCE OF. 
□ 20 BACTERIOPHAGE \ DNA 

' Lambda phage D N A is 4 3 ,502 base pairs long. It = 
sequence has been cczpieteiy determined, we have 
created one strand z: this as a :es: case in a compute 
simulation of the analysis. The table shows that z z e 
25 appropriate size of oligomer to use for a sequence :f 
this complexity is z h e 10-mer. With a matrix of ;I- 
mers, the size was 13 24 lines square. After 
"hybridisation" c : -.he lambda 10-mers in the co::::;:, 
46,277 cells were positive, 1957 had double' 
30 occurrences, 75 triple occurrences, and three quatrup 
occurrences. These ^5,577 positive ceils repress::; 
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^ known sequences, de:err.ined from their position in tne 

matrix. Each was expended by four x one base at tr.e 3' 
end and four x one : a s e at the 5 1 , end to give \ i x 
35 . 4 6,377 = 742,032 cells. This extended. set reduce: tne 
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number of double occurrences to 1 6 1 , a further 1 6 - 
fold extension brought the number down to 10, and one 
more provided a completely overlapped result. Of 
course, the same end result of a fully overlapped 
5 sequence could be achieved starting with a 4 1 ^ matrix, 
but the matrix would.be U c 0 0 times bigger than the 
matrix needed to represent all 10-mers, and most of the 
sequence represented on it would be redundant. 
5-2 LAYING DOWN THE MATRIX 
*Q The method described here envisages that the 

matrix will be produced by synthesising oligonucleotides 
in the cells of an array by laying down the precursors 
for the four bases, in a predetermined pattern, an 
example of which is described above. Automatic 
-5 equipment for applying the precursors has yet to be 
developed, but there are ,c':vious pos s'ibi 1 i t e s ; it 
should not be difficult t: adapt a pen plotter or other 
computer-controlled printing device to the purpose. 
The smaller the pixel size of the array the better, as 
20 complex genomes need very large numbers of cells. 

However, there are limits tc how small these can be 
made. 100 microns would be a fairly comfortable upper 
limit, but could probably not be achieved on paper for 
reasons of texture and diffusion. On a smooth 
-5 impermeable surface, sue r. is glass, it may be possible 
to achieve a resolution o f around 10 microns, for 
example by using a laser typesetter to preform a 
solvent repeilant grid, an: building the oligonucleo- 
tides in the exposed regit.-. s . One attractive 
- Dossibiiity, which allows abaooaticn of oresent 

techniques of oligc nucleotide s y n thesis, is to sinter 
micro porous glass in "i:r:!:cpic patches onto tne 
surface of a glass plate. Laying down very large 
number of lines or cots : : : 1 : take a long time, if tne 
- ^ printing mechanism were slew. However, a low cost ink- 



jet printer can print at speeds of about 10,000 spots 
per second* With this sort of speed/ 10 spots could 
be printed in about three hours. 

5 5.3 OLIGONUCLEOTIDE SYNTHESIS 

There are several methods of synthesising 
oligonucleotides. Most methods in current use attach 
the nucleotides to a solid- support of controlled pore 
size glass (CPG) and. are suitable for adaptation to 

10 synthesis on a glass surface. Although we know of no 
description of the direct use of oligonucleotides as 
hybridisation probes while still attached to the matrix 
on which they were synthesised, there are reports of 
the use of oligonucleotides as hybridisation probes on 

^5 solid supports to which they were attached after 

synthesis. PCT Application WO 85/01051 describes a 
jm^thod for synthesising oligonucleotides tethered to a 
CPG column. In an experiment performed by us,. CPG was 
used as the support in an Applied Bio-sytems 

2Q oligonucleotide synthesiser to synthesise a 13-mer 
complementary to the left hand cos site of phage 
lambda. The coupling steps were all close to 
theoretical yield. The first base was stably attached 
to the support medium through all the synthesis and 

25 deprotection steps by a covalent linX, 

5,4 ANALYSING SEVERAL SEQUENCES SIMULTANEOUSLY 

The method of this invention can be used to 
analyse several polynucleotide sequences 

2q simultaneously* To achieve this, the oligonucleotides 
may be attached to the support in the form of * (for 
example) horizontal stripes. A technique for doing 
this is described in Example 3 below. Each DISfA sample 
to be analysed is labelled and applied to the surface 

25 carrying the oligonucleotides in the form of a stripe 
(e.g. vertical) orthogonal to the oligonucleotide 
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stripes of the array. Hybridisation is seen at the 
intersections between oligonucleotide stripes and 

between them. 

5 Where sequence variations are known, an 

advantage of using this technique is that many 
. different mutations can be probed simultaneously by 
laying down stripes corresponding to each allelic 
variant. With a density of one oligonucleotide per 
■jq mm, and one "individual" per 5 mm, it should be 

possible to analyse 2000 loci on a plate *100 mm square. 
Such a high density of information where the 
oligonucleotides do identify specific alleles!, is not 
available by other techniques. 

15 

6. PROBES, HYBRIDISATION AMD DETECTION 

The yield of oligonucleotides synthesised on 

microporous glass is about 30 iimol/g. A patch of this 

material 1 micron thick by' 10 microns square would hold 
-12 ■ 

2 Q * 3 x 10 pmol', equivalent to about 2 g of human 
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D N A . The hybridisation reaction could therefore be 
carried out with a very large excess cf : r. e bound 
oligonucleotides over that in the probe. So it should 
be possible to design a s.ystem capable of 
5 distinguishing between hybridisation involving single 
and multiple occurrances of the probe sequence, as 
yield will be proportional to concentration at all 
stages in- t-h.e^ reaction. 

The polynucleotide sequence to be analysed may be 

10 of DNA or R N A . To prepare the probe, the polynucleotide 
may be' degraded to form fragments. Preferably it is 
degraded by a method which is as randec as possible, to 
an average length around the chosen lengtr. s of the 
oligonucleotides on the support, and oligomers of exact 

15 lengths selected by electrophoresis on a sequencing 

gel. The probe is then labelled. For example, 

oligonucleotides of length s may be end labelled. If 
32 

labelled-wit h" -P , the radioactive yield cf any 

individual s-mer even from total human DNA could be 
4 

20 more than 1C dpm/mg of total DMA. For detection, only 
a small fraction of this is needed in a patch ) 0-100 
microns square. This allows hybridisati:.-. conditions 
to be chosen to be close to the T m of duplexes, which 
decreases the yield of hybrid and decreases the rate cf 

25 formation, but increases the discriminating power. 

Since the bound oligonucleotide is in excess, signal 
need not be a problem even working close 
equilibrium. 

Hybridisation conditions can be c h c s e n to be those 
30 known to be suitable in standard procedures used z c 
hybridise zo filters, but establishing :::icu" 
conditions is important. In particular, lemperacure 
needs to be controlled closely, preferably :o better 
than +G.5"C. Particularly when the c noser, length of 
35 the oligonucleotide is small, the analysis needs to be 



• * 

-15- 

able to distinguish between slight differences of rate 
and/or extent of hybridisation. The equipment may need 
to be programmed for differences in base composition 
between different oligonucleotides. In constructing 
5 the array, it may be preferable to partition this into 
sub-matrices with similar base compositions. ' This may 
make it easier to define the Tm which may differ 
slightly according to the base composition. 
The choice of hybridisation solvent is 

10 significant. When 1M NaCl is used, G:C base pairs are 
more stable than AiT base pairs. Double stranded 
oligonucleotides with a high G+Ccontent have a higher 
Tm than corresponding oligonucleotides with a high A+T 
content. This discrepancy can be compensated in 
various ways: the amount of oligonucleotide .laid down 
on the surface of the support can be varied depending 
on its nucleotide composition; ' or the computer used to 
analyse the data can be programmed to compensate for 
variations in nucleotide composition.: A preferred 
• 20 method, which can be used either instead of or in 
addition to those already mentioned, is to use a 
chaotropic hybridisation solvent, for example: a • 
quarternary or tertiary amine as 'mentioned above. 
Tetramethylammoniumchloride (TMAC1) has proved 

2 cj particularly suitable, at concentrations in the range 2 
M to 5.5 H, At TMAC1 concentrations around 3,-5 M to 4 
M, the T^ dependence on nucleotide composition is 
greatly reduced. 

The nature of the hybridisation salt used also has 

3q a major effect on the overall hybridisation yield. 

. Thus, the use of TKAC1 at concentrations up to 5 M can 
increase the overall hybridisation yield by a factor of 
30 or more (the exact figure depending to some extent 
on nucleotide composition) in comparison with 
hybridisation using 1M NaCl. Manifestly, this has 
important implications; for example the amount of 
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probe material that needs to be used to achieve a given 
signal can be much lower. 

Autoradiography, especially with 32 P causes image 
degradation which may be a limiting factor 'determining 
^ resolution; 'the limit for silver halide films is around 
25 microns. Obviously some direct detection system . 
would be better. Fluorescent. probes are envisaged; 
given the high concentration of the target 
oligonucleotides, the low sensitivity of fluorescence 

10 may not be a problem. 

We have considerable experience of scanning auto- 
radiographic images with a digitising scanner. Our 
present design is capable of resolution down to 25 
microns, which could readily be extended down to less 

15 than present application, depending on the quality of- 
the hybridisation reaction, and how good it is at 
distinguishing absence of a sequence from the presence 
of one or more. Devices for measuring astronomical 
plates have an accuracy around 1 jf. Scan speeds are 

2Q such that a matrix of several million cells can be 

scanned in a few minutes. Software for .the analysis of 
the data is straight-forward, though the large data 
sets need a fast computer- 
Experiments presented below demonstrate the 

25 feasibility of the claims. 

Commercially available microscope slides (BDH. :/ 
Super Premium 76 x 26 x 1 mm) were used as supports. 
These were derivatised with a long aliphatic linker 
that can withstand the conditions used for the 

30 
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deprotection of the aromatic heterocyc 1 i c bases , i.e. 

30% NH at 55° for 10 hours. The linker, bearing a 
3 

hydroxyl group which serves as a starting point for the 
subsequent oligonucleotide, is synthesised in two 

D steps. The slides are first treated with a 25% 

solution of 3-glycidoxypropyitriethoxysilane in xylene 
containing several drops of Hunig's base as a catalyst. 
The reaction is carried out in a staining jar, fitted 
with a drying tube, for 20 hours at 90 C. The slides 

10 are washed with MeOH, Et^O and air dried. Then neat 
hexaethy l'ene glycol and a trace amount of cone, 
sulphuric acid are added and the mixture kept at 80 
for 20 hours. The slides are washed with MeOH, Et 0 , 
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air dried and stored desiccated at -20 until use. 
15 This preparative technique is described in British 

Patent Application 8822228.6 filed 21 September 1988. 

The oligonucleotide synthesis cycle is performed 
as follows: 

The coupling sol ut i on" i s- made up fresh for each 

20 step by-mixing 6 vol. of C.5M tetrazole in anhydrous 

acetonitrile with 5 vol. cf a 0.2M solution of the 

required bet a- cy ancet hy 1 pho s pho r am i d i te . Coupling tim 

is three minutes. Oxidation with a 0 . 1 M solution 

" of I in THF/pyridine/H 0 yields a stable phospho- 
- 2 2 

20 ■ triester bond. Detri ty 1 = - ion of the 5' end . with 2% 

trichloroacetic acid in d i ch 1 oromethane allows further 

extension cf the oligonucleotide chain. There was nc 

capping step since the excess cf phos phoramid i tes usee 

over reactive sites on the slide was iarge enough tc 

30 drive the coupling to ccr.: 1st ion . Af-ter the synthesi 

is completed, the oligonucleotide is deprotected in 30 

NH for 10 hours at 55°. The chemicals used in the 
3 ■ *• ■ 

coupling step are moisture-sensitive, 2nd this critics 

step must be performed uncer anhydrous conditions in a 

55 sealed container, as follows. The shape of the patch 



10 



-18- 

to be synthesised was cut out of a sheet of silicone 
rubber (76 x 26 x 0.5 mm) which was sandwiched between 
a microscope slide, derivatised as described above, and 
a piece of teflcn of the same size and thickness. To 
this was fitted a short piece of plastic tubi'ng that 
allowed us to inject and withdraw the coupling solution 
by syringe and to flush the cavity with Argon. The 
whole assembly was held together by fold-back paper 
clips. After coupling the set-up was disassembled and 
the slide put through the subsequent chemical reactions 
(oxidation with iodine, and de.tr i ty 1 at ion by treatment 
with TCA) by dipping it into staining jars. 
EXAMPLE 1 . 

As a first example we synthesised the sequences 

cliqo-dT -oliqo-dT on a slide by qraduallv decreasina 
10 * 14 

the level of the coupling solution in steps 10 to 14. 

Thus the 10-mer was synthesised cn the upper part of 

the slide, the 14-mer at the bottom - and" the 11, 12 and 

13-mers were in between. We used 10 pmol oliao-dA 

■ 32 "12 
labelled at the 5 5 end with P by the polynucleotide 

kinase reaction to a total activity of 1.5 million 

c.p.m., as a hybridisation probe. Hybridisation was 

carried out in a perspex (Plexiclas) container made to 

-it a microscope slide, filled with 1.2 ml of 1M NaCl 

o 

in TE, 0.1% SDS, for 5 minutes at 20 . After a... short 
rinse in the same solution withi-Jt oligonucleotide, we 
were able to detect more than 2000 c.p.s. with a 
radiation monitor. An autorac :: graph. ..showed that all 
the counts came from .the area were the oligonucleotide 
■-ad been synthesised, i.e. there was no non-specific 
winding to the glass or to the region that had been 
cerivatised with the linker only. After partial 
elution in 0.1 M NaCl differential binding to the 
target is detectable, i.e. less binding tc the shorter 
3 5 t.nan the longer oli'go-dT. By gradually heating the 
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slide in the wash solution we determined the T (mid- 

om 

point c r transition when 50% eluted) tc be 33 . There 

o 

were no counts detectable after incubation at 39 . 
The hybridisation, and melting was repeated eight times 
^ with no diminution of the signal. The result is 
reproducible. We estimate that at least 5% of the 
input counts were taken up- by the slide at each cycle. 
EXAMPLE 2. 

In order to determine whether we would be able to 
10 distinguish between matched and mismatched oligo- 
nucleotides we syrvthesised two sequences 3' CCC GCC GCT 
GGA (cosL) and 3' CCC GCC TCT GGA, which' differ by one 
base at position 7. All bases'except the seventh were 
added in a rectangular patch. At the seventh base., 
' 5 half of the rectangle was exposed in. turn to add the 

two different bases, in two stripes. Hybridisation of 

cosR orcbe o 1 i aonuc 1 eot i de (5* GGG CGG CGA. CCT) (kinase 

32 - 
labelled with. P to 1.1. million c.p.m., 0.1 M NaCl, 

TE, C.11 SOS) was for 5 hours at 32 . The front of the 

20 slide showed 100 c.p.s. after rinsing. Autoradiography 
showed that annealing occurred only t: c h e part of the 
slide with the fully como 1 ement ary o 1 i z :n uc 1 eot i de . 
No signal was detectable on the patch with the 
mi smatcred sequence . 

2 5 EXAMPLE 3. 

For 2 further study of the effects of mismatches 
or shorter sequences on hybridisation behaviour, we 
constructed two arrays; one (a) of 2^ ,c \ : gonuc 1 eot i de s 
and the :ther (b) of 72 oligonucleotides. 

3° These arrays were set cut as shewn in Table" 1(a) 

and 1(b*, . The masks used to lay down these arrays were 
d i f f e r e r. *: from those used in previous experiments. 
Lengths :f silicone rubber tub i no ( Iran c.d.) were glued 
with silicone rubber cement to the surface of plain 

35 microsccce slides, in the form of a "U". Clamping 
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these masks against a derivatised microscope slide 
produced a cavity into which the coupling solution was 
introduced through a syringe. In this way cnly the 
part of the slide within the cavity came into contact 
with the phosphor am idite solution. Except in the 
positons of the mismatched bases, the arrays listed in 
Table 1 were laid down using a mask which covered most 
of the width of the slide. Off-setting this r.ask by 
3mm up or down the derivatised slide in subsequent 
coupling reactions produced the o 1 i gnuc 1 eot i des 
truncated at the 3' or 5' ends. 

For the introduction of mismatches a -mask was used 
which covered half (for array (-a)) or one third (for 
array (b)) of the. width of the first mask. Th e bases 
^5 at positions six and seven were laid down in two of 

three longitudinal stripes. This led. to the synthesis 
of oligonucleotides differing by one base on each half 
(array (a)) or third (array (b)) of the slide. In 
other positions, the .sequences differed from the 
longest sequence by the absence of bases at the ends. 

In array (b), there were two columns of sequences 
between those shown in Table 1(b), in which the sixth 
and seventh bases were missing in all positions, 
because the slide was masked in a stripe by t ~ e 
silicone rubber seal. Thus there were a total cf.72 
different sequences represented on the slide i" 90' 
different positions. 

The 19-mer 5' CTC CTG AGG AGA AGT CTG C *=s used 
for hybridisation (2 million cprr, , 1 .2 'ml 0.'^ \ = C1 in 
30 TE, 0.1% SDS, 20 C ) . 

The washing and elution steps were relieved by 
autoradiography. The slide was kept in the * = 3 h i n g 
solution for 5 r.in at each elution step and trer- 
exposed (45 rr.ir. , intensified). Elution ter.pe^atures 

■3" .o 

were 23, 36, 4^, 47-,-55 -and 60 C respectively. 
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As indicated in the table, the oligonucleotides 

showed different meltina behaviour. Short ol i go- 

^o 

nucleotides melted before longer ones, and at 55 C, 
only the perfectly matched 19-mer was stable', all oth 
oligonucleotides had been eluted. Thus the method c 
differentiate between a 18-mer and a 19-mer which 
differ only by the absence of one base at the end. 
Mismatches at the end of the oligonucleotides and at 
internal sites can all be melted under conditions whe 
the perfect duplex remains. 

Thus we are able to use very stringent 
hybridisation conditions that el iminate .-annealing to 
mismatch sequences or to oligonucleotides differing i 
length by as little as one base. No other method usi 
hybridisation of oligonucleotides bound to the solid 
supports is so sensitive to the effects of mismatchin 

EXAMPLE 4. 

To test the application of the invention to 
diagnosis of inherited diseases, we hybridised the 
array (a), which carries the oligonucleotide sequence 
specific for the wild type and the sickle cell mutavl 
of the P -globin gene, with a 110 base pair fragment 
ON A amplified frcm-the p -clobin gene by means of the 
polymerase chain reaction (PCR). Total DNA from th* 
blood of a normal individual (1 microgram) was 
amplified by PCR in the presence of appropriate pri.m^ 
oligonucleotides.- The resulting 110 base pair frag- 
was purified by elecrcchores i s through an agarose ge 
After elution, a small sample (ca. 10'picogram) was 
labelled by using -""P-dCTP (50 microCurie) in a 
second PCR reaction. This PCR contained only the 
upstream priming oligonucleotide. After 60 cycles c 
amplification with an extension tine of 9 m i n . the 
product was removed fr:r: precursors by gel filtratic 
Gel electrophores i s'* c- -the radioactive product showe 
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major band corresponding in length to the 110 base 

fragment. One quarter of this oroduct (100,000 c.p.m. 

in 0.9 M NaCl, TE, 0.1% SDS) was hybridised to the 

o 

array (a). After 2 hours at 30 ca. 15000 c.p.m. had 
been taken uo. The melting behaviour of the hybrids 
was followed as described for the 19-mer in example 3, 
and it was found that the melting behaviour was similar 
to that of the oligonucleotide. That is to say, the 
mismatches considerably reduced the melting temperature 
of the hybrids, and conditions were readily found such 
that the perfectly matched duolex rema i ned -wherea s the 
mismatched duplexes had fully relted. 

Thus the invention can be used- to analyse long 
fragments of DNA as well oligonucleotides, and this 
example shows how it may be used to test nucleic acid 
sequences for mutations. In particular it shows how it 
may be applied to the diagnosis of genetic diseases. 
EXAMPLE 5. 

To test an automated system for laying down the 

20 precursors, the cosL oligonucleotide was synthesised 
with 11 of the 12 bases added in the way described 
above. For the addition of tre seventh base, however, 
the slide was transferred into an Argon filled chamber 
containing a pen plotter. The pen of the plotter had 

25 been replaced by a component, -"abricated from Nylon, 

which had the same shape and dimensions as the pen, but 
which carried.a pc lytetraf 1 ucr: ethyl ene (PTFE) tube, 
through which chemicals could :e delivered to the 
surface of the glass slice wh::h lay on the bed of the 

30 plotter. A microcomputer was 'jsed to control the 
plotter and the syringe pump *Mch delivered the 
chemicals. The pen, carrying zhe delivery tube from 
the syringe, was moved into : c e : t i o n above the slide, 
t h e pen was lowered and the p-~o activated to lay down 

35 ecu Dlina solution. ■*■ Filling -~e pen successively with 
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G, T and A phosphoramidite solutions an array of twelve 
spots was laid down in three groups of four, with three 
different oligonucleotide sequences. After 
hybridisation to cosR, as described in Example 2, and 
autoradiography signal was seen only over the four 
spots of perfectly matched oligonucleotides, where the 
dG had been added. 
E XAMPLE 6 

This example demonstrates the technique of 
analysing several DNA sequences simultaneously. Using 
the technique described in Example 3, a slide ; was 
prepared bearing six parallel rows of oligonucleotides 
running along its length. These comprised duplicate 
hexadecamer sequences corresponding to antisense 
1^ sequences of the p-globin wild-type (A) , sickle ceil - 
(S) and C mutations. 

Clinical samples of AC, AS and SS DNA were 
procured. Three different single -stranded probes of 
110 nt length with, approx. 70,000 c.p.m, in 100 p.1 1K 
2Q -NaCl", TE pH 7.5, 0.1% SDS, viz AC, AS y and SS; DNA were 
prepared. Radiolabeled nucleotide was included in the 
standard PGR step yielding a double-stranded labelled 
fragment. It was made single-stranded with 
Bacteriophage \ exohuclease that allowed to selectively 
2g digest one strand bearing a 5' phosphate. This was 

made possible by phosphoryiating the downstream primer 
with T4 Polynucleotide kinase and ('cold') ATP prior to 
PGR. These three probes were applied as three stripes 
orthogonal to the surface carrying the six 
oligonucleotide stripes. Incubation was at 30 *C for 2 
hours in a moist chamber. The slide was then rinsed at 
ambient temperature, then 45 "C for 5 minutes and 
exposed for 4 days with intensification. The genotype 
of each clinical sample was readily determined from the 
autoradiographic signals at the points of intersection. 

35 
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EX AMF EE 7 . • 

A plate was prepared whose surface carried an 
array of all 256 octapurines. That Is to say| # the 
array comprised 256 oligonucleotides each consisting 
5 of a different sequence of A and G nucleotides. This 
array was probed with a mixture comprising all 256 
octapyrimidines, each end labelled by means of 
polynucleotide kinase and y- P-ATP. Hybridisation 
was performed for 6-8 hours at 4*C, 

10 * n consecutive experiments the hybridisation 

solvent was changed through the series 1M NaCl 
(containing lOmM Tris.HCi pH 7-5', 1mM EDTA, 7j% 
sarcosine) and 2M, 2.5M, 3M, 3.5M, 4M, 4.5M, 5K and 
5.5M TMAC1 (all containing 50mM Tris.HCi pH 8..0, 2mM 

^ EDTA, SDS at less than 0*04 mg/ml) . The plate was 

rinsed for 10 minutes at 4*C in the respective solvent 
to remove only loosely matched molecules, sealed in a 
plastic bag and exposed to a Phorphorlmager storage 
phosphor screen at 4*C overnight in the dark. 

2Q The following table quotes relative signal 

intensities, at a given salt concentration, of hybrids 
formed with oligonucleotides of varying a content* In 
this table, the first row refers to the oligonucleotide 
GGGGGGGG, and the last row to the oligonucleotide 

2^ AAAAAAAA. It can be seen that the difference in 

response of these two oligonucleotides is marked in 1M 
NaCl, but much less marked in 3M or 4M THAC1 . 
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Relative Intensities at give n Salt: Concentration 



Solvent 



1M NaCl 
2M TMAC1 
3M TMAC1 
4M TMAC1 



Number of A's 
0 4 



100 
100 
, 70 
60 



30 
70 
100 
100 



20 
30 
40 
40 



15 



20 



The following table indicates relative signal 
intensities obtained, with octamers containing 4A*s and 
4G's, at different hybridisation salt concentrations. 
It can be seen that the signal intensity is 
dramatically increased at higher concentrations of 
TMAC1 . 

Relative Intensities at di fferent Salt Concentrations 



25 



Solvent Yield of hybrid 



• 1M NaCl 100 

2M TMAC1 200 

3M TMAC1 700 

4M TMAC1 ' 2000 
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In conclusion, we have demonstrated the following: 
1 . It is possible to synthesise oligonucleotides in 
good yield on a flat glass plate. 
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2, Multiple sequences can be synthesised on! the sample 
in small spots, at high density, by a simple manual 
procedure, or automatically using a computer controlled 
device . 

g 3. Hybridisation to the oligonucleotides ori the plate 
can be carried out by a very simple procedure;* 
Hybridisation is efficient, and hybrids can be detected 
by a short autoradiographic exposure . 

4. Hybridisation is specific. There is no jdetectable 

, i 

1Q signal on areas of the plate where there. are !no 

oligonucleotides. We have tested the effects of mis- 
matched bases, and found that a single mismatched base 
at any position in oligonucleotides ranging in length 
from 12-mer to 19-mer reduces the stability df the 

15 hybrid sufficiently that the signal can be reduced to a 
very low level, while retaining significant 
hybridisation to the perfectly matched hybrid. 

5. The oligonucleotides are stably bound to: the glass 
and plates can be used for hybridisation repeatedly. 
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The invention thus provides a novel way of 



analysing nucleotide sequences, which should find a 
wide range of application. We list a number of 
potential applications below: 
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Small arrays of oligonucleotides as fingerprinting 

and mapping tools 

Analysis of known mutations including genetic diseases. 

Example 4 above shows how the invention may be 
5 used to analyse mutations. There are many 

applications for such a method, including the detection 
of inherited diseases. 

Genomic fingerprinting. 

In the same way as mutations which lead to disease 

10 can be detected, the method could be used to detect 

point mutations in any stretch of DNA. Secuences are 
now available for a number of regions containing the 
base differences which lead to rest r i c t i on fragment 
length polymorphisms (RFLPs). An array of oligo- 

'5 nucleotides representing such polymorphisms could be 
made from pairs of oligonucleotides representing the 
two allelic restriction sites. Amplification of the 
sequence containing the RFLP , f o 1 1 owed by hybridisation 
to the plate, would show which alleles were present in 

2G the test genome. The number of oligonucleotides that 
could be analysed in a single analysis could be quite 
Large. Fifty pairs made from selected alleles would 
be enough t: give a fingerprint unique to an individual. 
Linkage ana lysis. 

?5 • Applying the method described in the last 

paragraph iz a pedigree would pinpoint reccr.: i nations'. 
Each pair of spots in the array would give the 
information -hat is seen in the track c* the RFLP 
analysis, using gel electrophoresis and' b 1 z z z i ng , that 

3° is now routinely used for linkage studies. It should 
be possible zo analyse many alleles in a single 
analysis, by hybridisation to an array of allelic pairs 
of oligonucleotides, greatly simplifying the methods 
used to fin: linkage between a DNA polymorphism and 

35 phenotypic -=rker such as a disease gene. 

The ex er. pies above could be carried-cut using the 
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method we have developed and confirmed by experiments. 
Large arrays of oligonucleotides as sequence 

reading tools. 

We have shown that oligonucleotides can be 
synthesised in small patches in precisely determined 
positions by one of two methods: by delivering the 
precursors through the pen of a pen-plotter, or by 
masking areas with silicone rubber. It is obvious how 
a pen plotter could be adapted to syntnesise large 
arrays with a different sequence in each position. Fcr 
some applications the array should be a predetermined, 
limited set; for other applications, the array should 
comprise every sequence of a predetermined length. The 
masking method can be used for the latter by laying. 
^ down the precursors in a mask which produces 

intersecting lines. There are many ways in which this 
can be done and we give one example for illustration: 
1. The first four' bases ,' A ,' C . G, T, are laid in r'our 
broad stripes on a square plate. 
20 2 . The second set is laid down in four stripes equal 
in width to the firs-., and orthogonal to them. The 
array is now composed of all sixteen di nuc leoti des . 

3. The third and "fourth layers are laid down in f:ur 
sets of four stripes one quarter the width of the first 

25 stripes. Each set cf four narrow stripes runs ..wit' in 
one of the broader stripes. The array is now composed 
of all 256 tetranuclectides . 

4. The process is. repeated, each time laying down two 
layers with stripes which are one quarter the width cf 

30 the previous two layers. Each layer added increases 

the length of the oligonucleotides by one base, an: the 
number cf d i f f er e.nt -c 1 i ccr.uc 1 eot i de sequences by a 

factor of four. 

The dimensions cf such arrays are determined c> 
35 the width of the stripes. The narrowest ...stri pe we 



have laid is 1mm, but this is clearly, not the lowest 
limit. 

There are useful applications for arrays in which 
part of the. sequence is predetermined and part made up 
of all possible sequences. For example: 

Characterising m R N A populations. 

Most mRNAs in higher eukaryotes- have the sequence 

AAUAAA close to the 3' end. The array used to analyse 

mRNAs would have this, sequence all over the plate. To 

analyse a mRNA population it would be hybridised to an 

array composed of all sequences of the type N AATAAAN . 
_ m n 

For m + n = 8, which should be enough to give a unique 

oligonucleotide address to most" of the several thousand 

mRNAs that is estimated to be present in a source such 

as a mammalian cell, the array would be .256 elements 

square. ■ The 256 x 255 elements wo,uld be laid on the 

AATAAA using the masking method described above. With 

stripes of around 1mm, the array would be ca. 256mm 

square.. 

This, analys is would measure the complexity of the 
mRNA population and could be used as a basis for 
comparing populations from different cell types. The 
advantage of th i s approach is that the differences in 
the hybridisation pattern would provide the sequence of 
oligonucleotides that could be used as probes to... 
isolate all the mRNAs that differed in the populations. 

Sequence determination. 

To extend the idea to determine unkown sequences, 
using an array compos.ed of all possible oligonucleotide 
of a chosen length, requires larger arrays thanwe have 
synthesised to date. However, it is possible to scale 
down the size of spot and scale up the numbers to those 
required by extending the -.ethods we have developed and 
tested on small arrays. Cur experience shows that the 
method is much simpler -in cperation than the_gel based 
methods. 



TABLE 1 



For 


Examples 


3 and 4 


array (a) 


was set cut as 


fol 


lows 








20 


GAG 


GAC 


TCC TCT 


ACG 




20 


GAG 


GAC 


aCC 


TCT ACG 






35 


GAG 


GAC 


TCC TCT 


GAC G 




20 


GAC 




aCC 


TCT 


GAC 


G 




36 


GAG 


GAC 


TCC TCT 


AGA CG 




20 


GAC 


GAC 


aCC 


TCT 


AGA 


CG 




47 


r~ « " 

U/nO 


GAC 


TCC TCT 


CAG ACG 




36 


GAG 




aCC 


TCT 


CAG 


ACG 




60 


GAG 


GAC 


TCC TCT 


TCA GAC 


G 


47 


GAG 


GAC 


aCC 


TCT 


TCA 


GAC 


G 




.AG 


GAC 


TCC TCT 


TCA GAC 


G 


42 


.AG 


GAC 


aCC 


TCT 


TCA 


GAC 


G 


56 


..G 


GAC 


TCC, TCT 


TCA GAC 


G 


42 


..G 


GAC 


aCC 


TCT 


TCA 


GAC 


G 


47 




GAC 


TCC TCT 


TCA GAC 


G 


42 




GAC 


aCC 


TCT 


TCA 


GAC 


G 


42 




.AC 


TCC TCT 


TCA GAC 


G- 


36 




.AC 


aCC 


TCT 


TCA 


GAC 


G 


36 




..C 


TCC TCT 


TCA GAC 


G 


36 




..C 


aCC 


TCT 


TCA 


GAC 


G 


35 






TCC TCT 


TCA GAC 


G 


36 






aCC 


TCT 


TCA 


GAC 


G 


36 






.CC TCT 


TCA GAC 


G 


36 






.CC 


TCT 


TCA 


GAC 


G 



For 
2G vjAj 
20 GAG 
20 GAG 
20 GAG 
20 GAG 
20 GAG 
20 GAG 
20 GAG 

32 gag 

22 GA2 
42 GAG 
52 GAG 
42 .AG 
42 ..G 
37 ... 
32 ... 
32 ... 
32 ... 

3 r l '* 

cciu 
eve" 
(a.t 



example 3 array (b) was set out as follows: 



GAt TC 
GAt TCC 
GAt TCC T 
GAt TCC TC 
GAt TCC TCT 
'GAt TCC TCT T 
GAt TCC TCT TC 
GAt TCC TCT TCA 
GAt TCC TCT TCA G 
GAt TCC TCT TCA GA 
GAt TCC TCT TCA GAC 



20 GAG GAC TC 
20 GAG GAC TCC 
20 GAG GAC TCC T 
20 GAG GAC TCC TC 
20 GAG GAC TCC TCT 
20 GAG GAC TCC TCT T 
20 GAG GAC TCC TCT TC 
20 GAG GAC TCC TCT TCA 
42 GAG GAC TCC TCT TCA G 
47 GAG GAC TCC TCT TCA OA 
52 CAG GAC TCC TCT TCA GAC 



20 GAG GAC aC 
20 GAG GAC aCC 
20 GAG GAC aCC T 
20 GAG GAC aCC TC 
20 GAG GAC aCC TCT 
20 GAG GAC aCC TCT T 
20 GAG GAC aCC TCT TC 
20 GAG GAC aCC TCT TCA 
20 GAG GAC aCC TCT TCA G 
32 GAG GAC aCC TCT TCA GA 
42 GAG GAC aCC TCT TCA GAC 



•GAt TCC TCT TCA GAC G 60 GAG GAC TCC TCT TCA. GAG G 52 GAG GAC aCC TCT TCA. GAC G 

GAt TCC TCT TCA GAC G 52 .AG GAC TCC TCT TCA GAC G 42 .AG GAC aCC TCT TCA GAC G 

GAt TCC TCT TCA GAC G 52 ..G GAC TCC TCT TCA. GAC G 42 ..G GAC aCC TCT TCA GAC G 

GAt TCC TCT TCA GAC G 47 ... GAC TCC TCT TCA GAC G' 37 

..At TCC TCT TCA GAC G -42 . ; . .AC TCC TCT TCA GAC G 32 

. .t TCC TCT TCA GAC G 42 C TCC TCT TCA GAC G'32 

... TCC TCT TCA GAC G 32 TCC TCT TCA GAC G 32 

e=n the three column? of array (b) lis-.ed above, were two 
-.is, in which bases 6 and 7 from the left were missing in 

line. These sequences all melted a - . 20 or 32 degrees. 
; mismatch base (.) missing base. . 



; .GAC aCC TCT TCA GAC G 

. .AC aCC TCT TCA GAC G 

. ..C aCC TCT TCA GAC G 

. ... aCC TCT TCA GAC G 



